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Abstract. We study high-energy proton-nucleus collisions within the Colour Glass 
Condensate framework, and compute the probabilities of having a definite number of 
scatterings in the nucleus with a momentum transfer larger than a given cut. Various 
properties of the distribution in the number of multiple scatterings are investigated, 
and we conclude that events with monojets are very unlikely, except for extreme values 
of the saturation scale Q s . 



1. Introduction 

High-energy collisions of hadrons probe partons with very small x, whose occupation 
number becomes large, giving rise to the phenomenon of parton saturation [1]. The 
saturation is enhanced in a large nucleus, in which the nucleon wave functions overlap, 
and is expected to be revealed in high-energy collisions involving such a heavy nucleus. 

Different frameworks implementing parton saturation have been developed, 
especially the McLerran-Venugopalan (MV) model [2,3] and later the so-called "Colour 
Glass Condensate" (CGC) [4]. A general feature in these approaches is the description 
of small- x partons as classical Yang-Mills fields, radiated by the static large- x colour 
sources. Computations of particle production in the collision of two saturated nuclei 
are then rather involved [5] and require numerical approaches. However, in interactions 
involving one small projectile, the colour sources describing the latter can be treated as 
weak: particle production can then be calculated by considering the relevant amplitude 
only at lowest order in this source, which allows one to obtain analytical expressions [6]. 

Hereafter, we shall focus on single-gluon production [7]. While it is known that 
the CGC framework automatically includes all multiple scatterings that lead to the 
production of a gluon with a high transverse momentum k±, we shall study in further 
detail the distribution in the number of these scatterings [8], and more specifically, the 
scatterings with a momentum transfer larger than an arbitrary threshold fc™ m . 



2. Models for the proton— nucleus interaction 



The number of gluons produced per unit of transverse momentum and per unit of 
rapidity in the collision between a proton and a saturated nucleus [7] involves the 
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Fourier transform C(kj_) of a correlator of Wilson lines [9]. When the distribution of 
colour sources inside the nucleus only has Gaussian correlations, as is the case in the two 
models we shall consider hereafter, this correlator can be written in a form that has an 
interpretation in terms of independent multiple scatterings a la Glauber. In particular, 
the contribution to the gluon yield of cases in which the outgoing gluon has undergone 
exactly n scatterings can be isolated. One can thus show that the conditional probability 
that a gluon that comes out with a momentum kj_ has undergone n scatterings with 
momentum transfers larger than a threshold /c™ in (and an arbitrary number of softer 
scatterings) reads [8] 
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where p is the number density of scattering centres in the nucleus, L is the longitudinal 
size of the nucleus, p% = pL is the density of scattering centres per unit of transverse 
area, and o"(kj_) is the differential cross-section of a gluon with a scattering centre (the 
integral of which is <7 to t)- We shall present results obtained with two versions of this 
cross-section: either that, proportional to kj 4 , of the MV model — in which the colour 
sources have a local Gaussian distribution — , or the slightly more complicated cross- 
section, oc (Ql/kj_) ln[l + {Ql/kY) 1 } with 7 0.64, of the non-local Gaussian effective 
theory ( "asymptotic model" ) that describes the gluonic content of a nucleus evolved to 
very small x values [10]. These two models also yield different relationships between the 
saturation scale Q s and the scattering-centre density pi in equation (1). 

3. Results 

Instead of computing directly the probability P n , equation (1), for each successive integer 
n, we rather calculate the generating function 

+00 

F(z,k ± \kf n ) = Y, P n( k i-\ k f n )z n , (2) 

n=0 

from which it is easy to extract the individual probabilities. This function, which has a 
quite simple analytic expression (equation (8) in reference [8]), also gives easy access to 
the average number of scatterings above the threshold fc™ n , defined as 

+00 

N(k ± \kf a ) =J2 nP n(k±\kT n ), (3) 

n=l 

which is obviously given by the value at z — 1 of the first derivative of equation (2). 
The distribution of the probabilities P n (k±\k± in ) and the average number of recoils 
N(k±\k™ m ), computed within the MV model, are displayed in figure 1. First, the width 
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Figure 1. Left: distribution of the probabilities to have n recoils when the produced 
particle acquires k± = 10 GeV, for various values of the threshold fc™ ln . Right: number 
of recoils vs. the threshold fc™ in , for various values of the acquired momentum k±. 
Both panels assume the MV model with Q 2 S =2 GeV 2 . 

of the multiplicity distribution decreases with increasing fc™ m , which is quite natural, 
since it is less and less likely to have events in with a large number of recoils when fc™ in 
increases. One also sees that for all values of the threshold such that Q s -C fc™ m < k±, 
the most likely number of recoils is n — 1, while n = is most likely when k± < fc™ in . 
This is paralleled by the fact that the number of recoils is very close to one for any value 
of such that Q s -C fc™ m < k±: when the gluon acquires a large momentum k±, there 
is always one, and only one, hard recoil that provides most of this momentum. On the 
other hand, iV(/c_i_|/c™ m ) grows significantly at small fc™ m , where it does not depend on 
k±. This behaviour can be understood, since the generating function can be computed 
analytically when k± is much larger than the other two scales Q s and fc™ m . One then 
finds that the probabilities P n follow a Poisson distribution shifted by one, corresponding 
to the compulsory hard scattering that provides most of the acquired momentum; the 
average value of this distribution is controlled by Q s and fc™ in only, independent of k±. 
This analytic prediction is represented by the solid line in the figures. 

The interpretation in terms of a single hard scattering that provides most of the 
momentum is further supported by the momentum distribution of the recoils, left panel 
of figure 2. Indeed, this distribution consists of a universal semi-hard component made 
of recoils with momenta < Q s , and a second component peaked around fc™ m = k±. In 
other words, the most likely configuration is that of a di-jet, the high-momentum gluon 
with k± being balanced by a second hard parton with similar momentum, rather than 
a monojet, which would correspond to having many soft partons on the away side. 

When the saturation scale Q s is varied (figure 2, right panel), the second, hard 
component around /c™ in = k± is unchanged — as it should, since it is controlled by 
hard physics, not by saturation. On the other hand, the semi-hard region is affected 
by changes in Q s , with more and more scatterings as Q s increases. This growth in the 
number of soft recoils, which is quite spectacular within the MV model (open symbols 
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Figure 2. Distribution of the number of recoiling centres vs. the threshold fc™ m . 
Left: for various values of the acquired momentum kj_, within the MV model with 
Q 2 S = 2 GeV 2 . Right: dependence on the saturation Q s within the MV (open symbols) 
and asymptotic (full symbols) models, for a fixed value k± = 10 GeV. 



in figure 2), is however significantly reduced in the asymptotic model (full symbols). 
This lesser sensitivity is believed to reflect the presence of leading-twist shadowing in 
the latter model: first, the shadowing somehow "hides" the scattering centres from the 
passing gluon. In addition, since the effect of piling up more and more colour sources 
in the nucleus is tamed by the shadowing, the dependence on Q s is weaker in the 
asymptotic model than in the shadowing-free MV model. 

Note also that the two-component structure of the momentum distribution of the 
recoils is less marked in the asymptotic model than in the MV model (figure 2, right 
panel). Especially, the dip between the two regions is filled up. This means that one 
should expect that the distribution of momenta in the "away-side jet" will become 
flatter as one probes the nucleus at smaller and smaller x values. The di-jet structure, 
however, still persists, unless the value of the saturation momentum Q s becomes very 
large, opening the possibility of having a hard jet balanced by a large number of softer 
particles on the other side. 
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